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ABSTRACT 

We present the first phase in our ongoing work to use Sloan Digital Sky Survey 
(SDSS) data to create separate white dwarf (WD) luminosity functions for two or 
more different mass ranges. In this paper, we determine the completeness of the 
SDSS spectroscopic white dwarf sample by comparing a proper-motion selected 
sample of WDs from SDSS imaging data with a large catalog of spectroscopically 
determined WDs. We derive a selection probability as a function of a single color 
(g — i) and apparent magnitude (g) that covers the range —1.0 < g — i < 0.2 and 
15 < g < 19.5. We address the observed upturn in log g for white dwarfs with 
Teff < 12,000K and offer arguments that the problem is limited to the line profiles 
and is not present in the continuum. We offer an empirical method of removing 
the upturn, recovering a reasonable mass function for white dwarfs with T^g < 
12,000K. Finally, we present a white dwarf luminosity function with nearly an 
order of magnitude (3,358) more spectroscopically confirmed white dwarfs than 
any previous work. 
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function 



^The University of Texas at Austin, Department of Astronomy, 1 University Station C1400, Austin, TX 
78712-0259 

^Instituto de Fisica, Universidade Federal do Rio Grande do Sul, 91501-900 Porto- Alegre,RS, Brazil 
^Gemini Observatory, Hilo, HI 96720, USA 

^Institut fiir Theoretische Physik und Astrophysik, Universitat Kiel, 24098 Kiel, Germany 

^Facultad de Ciencias Astronomicas y Geoffsicas, Universidad Nacional de La Plata, Paseo del Bosque 
S/N, 1900, La Plata, Argentina 

^Instituto de Astrofisica La Plata, lALP, CONICET 



- 2 - 



Introduction 



Because white dwarfs cannot replenish the energy they radiate away — any residual nu- 
clear burning is negligible and gravitational contraction is severely impeded by electron 
degeneracy — their luminosity decreases monotonically with time. A thorough knowledge 
of the rate at which WDs cool can provide a valuable "cosmic clock" to determine the 
ages of many Galactic populations, includ ing the disk (IWinget et al.l 119871 : iLiebert et al. 

Knox et aDll999h. and open and globular clusters Jciaveil 



1988 



Leggett et al. 



1998 



von 



2006 



Hansen et al. 



2007 



1995 



lippel et al.lll995l:lRicher et al.lll998l:lciaver et all boOll: lllansen et al1l2002l . 120041 : Ivon Hippel et al 



Jeffery et al.l 120071 ). With more accurate models of the cooling 



physics of white dwarfs, heavily constrained by empirical evidence, it may be possible to 
determine absolute ages with greater precision than using main-sequence evolution theory. 
In addition to applications in astronomy, white dwarfs allow us to probe the physics of 
degenerate matter at temperatures and densities no terrestrial laboratory can duplicate. 



as 



Attern p ts at an er npirical luminosity function (LF) for white dwarfs date as far back 
LuvtenI (119581) a nd IWeidemanru (119671). The low- lum inosity shortfall, discovered by 



Liebert et al.l (jl979l ). and attributed by lWinget et al.l (jl987l ) to the finite age of the Galactic 
disk, was confirmed and explored m ore fully when a gr e ater volume of reliable data on low- 
luminosity WDs became available (ILiebert et al.lll988l : IWoodlll992l ). More recently, Sloan 
Digital Sky Survey (SDSS) photometric data have been used to provide a much more detailed 
luminosity functi on with more than an order of magnitude more white dwarfs than previ- 



ously attempted ([Harris et al.ll2006l) . as well a new LF of a large sample of spectroscopically 
confirmed WDs ( Hu et alT2007 ). However, to date no one has published a well-populated 
luminosity function that does not include a wide range of masses and spectral types. Thus, 
much of the important physics of white dwarf cooling remains buried in the data. 

Until recently, empirical WD luminosity functions, especially those derived from stars 
with spectra, have been hampered by a limited volume of reliable data. This has forced a 
trade-off between the number of stars included in a sample, and their homogeneity; either 
a broad range of temperatures, masses, and spectral types must be used, or else the sample 
population of stars would be so small as to render reliable conclusions difficult. Recently, the 
situation has changed dramatically. Data from SDSS DR4 have yielded nearly 10,000 white 



their effective temperatures and surface gravities (| 


Kleinman et al. 


2004 


2004; 


Eisenstein et al. 


2006; 


Hiigelmever et al. 


2006; 


Kepler et al. 


20071). 



In a companion paper to be published shortly, we intend to focus on how th e WD 
cooli n g rate changes with WD mass. Theoretical work has been done in this area (jWood 



1992 



Fontaine et al.l 1200 ll ). but to date, attempts at creating an empirical LF to explore 
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the effects of mass have rehed on hmited sample sizes (iLiebert et aLll2005l ). In order to 
further isolate the effect of mass, we have chosen to study only the DA white dwarfs — white 
dwarfs which show only lines of hydrogen in their spectra — which comprise ~86% of all white 
dwarfs. 

In addition to helping unlock the physics of white dwarfs, creating luminosity functions 
for several mass bins can also help to disentangle the effects of changes in cooling rates 
from changes in star formation rates. A burst or dip in star formation at a given instant in 
Galactic history should be recorded in all of the luminosity functions, regardless of mass, and 
could be confirmed by its position across the various mass bins. For example, a short burst 
of increased star formation would be seen as a bump in each luminosity function, occurring 
at cooler temperatures in the higher mass LF (these stars, with shorter MS lifetimes, have 
had longer to cool). On the other hand, features intrinsic to the cooling physics of the white 
dwarfs themselves should be seen in places that correspond with the underlying physics, 
which may be earlier, later, or nearly concurrent across mass b ins. These effects in clude 
neutrino coolin g, crystallization, th e onset of convective coupling (IFontaine et al.ll200ll ). and 
Debye cooling (lAlthaus et al.l 120071 ). 



The current paper lays the groundwork for this analysis. In Section [2l we introduce 
the data, examining the methods used to classify spectra and derive quantities of interest 
(dominant atmospheric element, Tes, and log g). We also address the observed upturn in 
log g for DAs below T^g ~ 12,000K. We present several lines of reasoning that the upturn is 
an artifact of the line fitting procedure, and propose an empirical method for correcting the 
problem. Section [3] outlines the methods used to construct the luminosity and mass function 
and determine error bars. 

In Section HI we present an analysis of the completeness of our data sample. We use 
a well-defi ned sample of pro per-motion selected, photometrically determined white dwarfs 
in SDSS (IHarris et al.l |2006| ) to determine our completeness and derive a correction as a 
function of g — i color and g magnitude. Finally, in Section O we present our best luminosity 
and mass functions for the entire DA spectroscopic sample and discuss the impact of both 
our empirical log g correction and our completeness correction. 



2. The Data 



Our white dwarf data comes mainly from lEisenstein et al.l (120061 ) . a catalog of spectro- 
scopically id entified white dw arfs from the Fourth Data Release (DR4) of the Sloan Digital 
Sky Survey (lYork et al.l l2000l ). The SDSS is a survey of ~8,000 square degrees of sky at 
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high Galactic latitudes. It is, first and foremost, a redshift survey of galaxies and quasars. 
Large "stripes" of sky are imaged in 5 bands (u,g,r,i,z) and objects are selected, on the 
basis of color and morphology, to be followed up with spectroscopy, accomplished by means 
of twin fiber-fed spectrographs, each with separate red and blue channels with a combined 
wavelength coverage of about 3800 to 9200A and a resolution of 1800. Objects are assigned 
fibers based on their priority in accomplishing SDSS science objectives, with high redshift 
galaxies, "bright red galaxies" and quasars receiving the highest priority. Stars are assigned 
fibers for spectrophotometric calibration, and other classes of objects are only assigned fibers 
that are left over on each p late. More detailed desc ripti ons of the targ e t sele ction and tiling 
algorithms can be found in IStoughton et al.l (120021 ) and iBlanton et al.l (120031 ) . 



Though white dwarfs are given their own (low priority) category in the spectroscopic 
selection algorithms, very few white dwarfs are targeted in this way. Rather, most of the 
white dwarfs in SDSS obtain spectra only th rough the "back door," most often when the 
imaging pipeline mistakes them for quasars. iKleinman et al.l (120041 ) list the various algo- 
rithms that target objects ultimately determined to be white dwarfs in DRl (their Table 1). 
White dwarfs are most commonly targeted by the QSO and SERENDIPITY_BLUE algo- 
rithms, with significant contributions also from HOT_STANDARD (standard stars targeted 
for spectrophotometric calibration) and SERENDIPITY_D1STANT. Of the significant con- 
tributors, the STAR_WHITE_DWARF category contributes the least to the population of 
WD spectra. 

The SDSS Data Release 4 contains nearly 850,000 spect ra. Several g r oups have al- 



ready attempted to sort through th e m to find white dwarfs: [Harris et al.l (120031 ) for the 



Early Data Release. iKleinman et al.l (120041 ) for Data Release 1 (DRl), and most recently. 



Eisenstein et al.l (120061 ) for the DR4, from which the majority of our data sample derives. 



Kleinman et al.l ( 


2004) 


from 


Eisenstein et al. 



(120061 ) with an expanded grid of models. A complete analysis of the 



methods by which candid ate objects are chosen, s pectra fitted, and quant ities of interest are 



calcul ated can be found in IKleinman et al.l (120041 ) , Eisenstein et al.l (120061 ) , and iKepler et al. 
( 120071 ). We put forth a brief outline here, with special attention paid to those aspects 
important to our own analysis. 

Objects in the SDSS spectroscopic database are put through s everal cuts in co l or de- 
signed to separate the WDs from the main stellar locus. Figure 1 in lEisenstein et al.l (120061 ) 
shows the location of these cuts. The chief failing of their particular choices of cuts, as 
noted by the authors, is that WDs with temperatures below ~8,000K begin to overlap in 
color-color space with the far more numerous A and F stars, and they have not attempted to 
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dig these stars out. The SDSS spectroscopic pipehne calculates a redshift for each object by 
looking for prominent lines in the spectrum. Objects with redshifts higher than z=0.003 are 
eliminated, unless the object has a proper motion from USNO-A greater than 0.3" per year. 
Because the spectroscopic pipeline is fully automated, occasionally DC white dwarfs show 
weak noise features that can be misinterpreted as low- confidence redshifts. Other types of 
WD, particularly magnetic WDs, can fool the pipeline as well. In the present paper we are 
concerned chiefly with DA white dwarfs, so this incompleteness is of importance only insofar 
as we use the entire set of white dwarf spectral types to derive our completeness correction, 
as outlined in Section |H We explore the implications of this more fully in that section. 



Eisenstein et al.l ( 120061 ) then use a minimization technique to fit the spectra and 



photometry o f the candidate ob j e cts with separate mo del atmospheres of pure hydrogen and 
pure helium (IFinley et al.l 119971 : iKoester et al.l 1200 ll ) to determine the dominant element, 
effective temperature, surface gravity, and associated errors. As their Figure 2 demonstrates, 
they recover a remarkably complete and uncontaminated sample of the candidate stars. They 
believe that they have recovered nearly all of the DA white dwarfs hotter than 10,000K with 
SDSS spectra. 

These stars form the core of our data sample. Their final table lists data on 1 0,088 white 
dwarfs. Of these, 7,755 are classified as single, non-magnetic DAs. iKepler et al.l (120071 ) re-fit 
the spectra for these stars using the same autofit method and Koester model atmospheres, 
but with a denser grid which also included models up to log g of 10.0. Where they differ from 
Eisenstein's, we use these newer fits in our analysis. Of these 7,755 entries, ~600 are actually 
duplicate spectra of the same star. For our analysis we take an average of the values derived 
from each individual spectra weighted by the quoted errors. Our final sample contains 7,128 
single, non-magnetic DA white dwarfs. 



As noted by iKleinman et al.l (120041 ) and others, the surface gravities determined from 
Sloan spectra show a suspicious upturn below temperatures of about 12,000K which increases 
at cooler temperatures, as shown in our Figure [H 

A number of separate pieces of evidence argue that this upturn in log g — and thus 
mass — is an artifact of the models and not a real effect. Not least among these is that no 
one has yet provided any satisfactory mechanism by which WDs could gain enough mass or 
shrink enough in radius as they cool to account for the magnitude of the effect. We do expect 
a slight increase in mass at cooler temperatures because in a galaxy of finite age, the cooler 
white dwarfs must come from higher mass progenitors. This is the reason for the upward 
slope of the blue dashed line in Figure [TJ However, this effect is clearly small compared to 
the upturn observed in the actual data. 
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Furthermore. lEngelbrecht &: Koesterl (120071 ). and lKepler et al.l (120071 ) demonstrated that 
the masses derived solely from the colors do not show an increase in mass for cooler stars, 
which indicates that the problem is not physical, but a result of either the line fitting pro- 
cedure or the line profiles themselves. 

Figure [3] further illustrates the above point. The upper panel sho ws the colors d erived 
from the synthetic spectra at the values of T^s and log g quoted by iKepler et al.l (120071 ) 
(i.e., the values in Figured]), overlaid on the actual SDSS photometry for the same objects. 
Contrast this with the lower panel, which instead shows the colors derived from the synthetic 
spectra when the excess log g has been removed (in a manner described below; the resulting 
values are shown in Figure [2]). The colors in the latter figure agree much better with the 
measured color of the object. 



Furthermore, iKepler et al.l (120071 ) found a similar increase in mean mass for the SDSS 
DB white dwarfs below Teg ~ 16,000K. They conclude that since a) the problem only shows 
up in the line profiles and not the continuum, and b) the onset of the effect in both hydrogen 
(DA) and helium (DB) atmosphere WDs occurs at just the effective temperature where the 
neutral species of the atmospheric element begins to dominate, then the problem lies in the 
treatment of line broadening by neutral particles. This is supported further by the fact that 
as the species continues to become more neutral (i.e., as the temperature drops), the problem 
grows worse. 

However, more recent model calculations indicate that neutral broadening is not impor- 
tant in the DA white dwarfs at temperatures down to at least than 8,500K. Other possible 
mechanisms to explain the observed upturn in log g include a flawed or incomplete treat- 
ment of convection, leading to errors in the temperature structure of the outer layers of 
the WD models, or the cony ective mixing of helium from a lower layer in the atmosphere 
( iBergeron et al.lll990l . Il995al ). The latter would require a hydrogen layer much thinner than 
any seismologically determined in a DA so far (iBradleyl Il998l . l200ll . |2006| ) . 



Until the problem with the model atmospheres is resolved, the best we can do is to 
empirically remove the log g upturn. For a given T^s, we subtract the excess in the measured 
mean value (as fit by the red solid lines in Figure [T]) over the theoretically expected mean 
(blue dashed line). Figure [2] shows the resulting values used. In fitting out the upturn this 
way, we make two implicit assumptions. First, that the excess log g is a function only of 
Tefj; if the problem is indeed due to the treatment of neutral particles, we would expect 
only a small dependence on log g. Second, we assume that the problem affects only the 
log g determination and not Teg. This latter assumption is unlikely to be true, as the two 
parameters are correlated. In Section [5] we explore more fully the impact of this fitting 
procedure on the luminosity and mass functions. 
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3. Constructing The Luminosity And Mass Functions 



Since we are dealing with a magnitude-limited sample, the most luminous stars in our 
sample can be seen to much further distances than the intrinsically fainter stars. We thus 
expect more of them, proportionally, than we would in a purely volume-lim ited sample, and 



must make a correction for th e different observing v olumes. As show n by IWood &: Oswalt 



Green 



mm andlGeiio et all (120061). the method of lSchmidtl fllQGSf ) (described more fully 



1980l : [Fleming et al.l Il986l ) provides an unbiased and reliable characterization 



m, e.g. 
of the WDLF. 

In the 1/Knax method, each star's contribution to the total space density is weighted in 
inverse proportion to the total volume over which it would still be included in the magnitude 
limited sample. Since the stars are not spherically distributed, but lie preferentially in the 
plane of the Galaxy, an additional correction for the scale height of the Galactic disk must 
be included. For the purposes of comparison with previous work, we adopt a scale height of 
250pc. 

To determine the absolu te magnitude of eac h WD, we use the effective temperatures 
and log g values provided by iKepler et al.l (120071 ) — as corrected in Section [2] — and fit each 
WD with an evolutionary model to determine the r nass and radius . For 7.0 < log g < 9.0, 
we use the mix e d C/0 models of IWoodI (Il995l ) and iFontaine et al.l (120011). as calculated by 
Bergeron et al.l (Il995bl ). For 9.0 < log g < 10.0, we use the models of lAlthaus et al.l (120051 ) 
with 0/Ne co res, including additional sequences for masses larger than 1.3 Mq calculated 
specifically for lKepler et al.l (120071 ). Once we know the radius, we can calculate the absolute 
magnitude in each Sloan band by con volving the synthetic WD atmospheres of Koester 
( iFinley et al.l Il997l: iKoester et al.ll200lf) with the Sloan filter curves. We apply bolometric 
corrections from iBergeron et al.l (Il995bl ) to determine the bolometric magnitude. For the 
handful of stars (~80-100) with log g values outside the range covered by Bergeron's tables, 
we use a simple linear extrapolation. 

We then determine photometric distances to each star from the observed SDSS g magni- 
tude. SDSS, being concerned mostly with extragalactic objects, reports the total inte rstellar 
absorption along each line of sight from the reddening maps of ISchlegel et al.l (119981 ). Since 
the objects in our sample lie within the Galaxy, and most of them w ithin a few hundre d 
parsecs, they are affected by only a fraction of this reddening. Following iHarris et al.l (120061 ) . 
we therefore assume: 1) that objects within lOOpc are not affected by reddening, 2) objects 
with Galactic height \z\ > 250pc are reddened by the full amount, and 3) that the reddening 
varies linearly between these two values. The distances and reddening are then fit itera- 
tively from the observed and calculated absolute g magnitude. In practice, the reddening 
correction makes very little difference to the final LF (typical Ag values range from 0.01 to 
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0.05). 

We calculate error bars on the luminosity function using a Monte-Carlo simulation, 
drawing random deviates in T^s, log g, and each band of photometry from gaussian dis- 
tributions centered around the measured value. The standard devi ations in T^ff an d log q 



we use for this scattering are 1.2 times the formal errors quoted in lEisenstein et al.l (120061 ) 
(their own analysis, based on repeated autofit measurements on duplicate spectra of the 
same stars, suggests that the formal errors derived by their method are ~20% too small). 
The photometry errors come directly from the SDSS database. After scattering the param- 
eters in this way, we recalculate the LF. We then add in quadrature the standard deviation 
of each LF bin after 200 iterations and the counting error for each bin (the errors for each 
individual star — taken to be of the order of the star's l/Vmax statistical weight — summed in 
quadrature). 

At a S/N of 16 — the mean for the stars in our sample brighter than g = 19.5 — formal 
errors in T^s and log g are of order 1.5%. When propagated through our code, the mean 
errors in Mhoi and mass are 0.35 dex and 9% (0.05 Mq) respectively. For the stars brighter 
than g = 19.0 used to compile our mass functions the average S/N is 19.5, leading to errors 
in Mboi and mass of 0.35 dex and 7% (0.04 Mq). 



4. Completeness Corrections 

The chief difficulty we have encountered in deriving our luminosity functions is un- 
raveling the complicated way in which SDSS objects are assigned spectral fibers. SDSS is 
foremost a survey of extragalactic objects and rarely targets white dwarfs for follow up spec- 
troscopy explicitly. Most of the objects in our sample are targeted by some other algorithm. 
In particular, there is considerable overlap in color between white dwarfs and many QSOs. 

A completeness correction could, in theory, be built from "first principles." We know, 
for each object in the SDSS spectroscopic database, by which algorithm(s) it was targeted (or 
rejected) for spectroscopy, and by which algorithm it was ultimately assigned a fiber. And 
for each algorithm, we know which objects were targeted, which were ultimately assigned 
a fiber, and which, of the targeted objects, turned out to be WDs. However, the selection 
process is a multi-variate function of 5 apparent magnitudes, and colors in spaces of as 
many as 4 dimensions (which vary based on the algorithm), as well as the complex tiling 
algorithm. We believe such an undertaking to be unnecessary. Instead we have ch osen to 



compare our sample with the stars used to derive the WDLF of iHarris et al.l (120061 ). Given 



certain assumptions about completeness and contamination in both data sets, we derive a 
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completeness correction as a function of a single color index [g — i) and g magnitude. 



The [Harris et al.l (120061 ) sample comes from photometric data in the SDSS Data Release 
3. They selected objects by using the reduced proper motion diagram to separate WDs from 
more luminous sub dwarfs of the same color. Briefly, they used color and proper motion 



(from USNO-B iMunn et al.l 120041 ) to determine WD candidates from SDSS imaging data. 
They then fit candidates with WD model atmosphere colors to determine temperatures and 
absolute magnitudes, from which they derived photometric distances and — together with 
proper motion — tangential velocities. In order to minimize contamination, they adopted a 
tangential velocity cutoff of 30 km/s and rejected all stars below this limit. The remaining 
6,000 objects are, with a high and well-defined degree of certainty (~ 98 — 99%), likely to 
be white dwarfs. 

If the database of SDSS spectra were complete, all of these objects would (eventually) 
have spectra, and all but the 1-2% of contaminating objects would be confirmed to be WDs. 
Furthermore, all of the WDs that did not make i t into the Harris et al. sample — because they 



were either missing from the iMunn et al.l (120041 ) proper motion catalog, or had a tangential 



velocity below 30 km/s — would also all have spectra. In such a perfect world, of course, 
no completeness correction would be necessary. However, since SDSS does not obtain a 
spectrum of every object in its photometric database, a significant percentage of the objects 
in Harris et al. will not have spectra, or else will be dropped at some later point by Eisenstein 
et al. and thus not make it into our spectroscopic sample. Our goal, then, is to look at all 
of the WDs in the Harris et al. sample that potentially could have made it into our sample, 
and determine which ones in fact did. If we assume that the WDs not in Harris et al. follow 
the same distribution (an assumption we discuss more fully below), then we can take this as 
a measure of the overall detection probability and invert it to get a completeness correction. 

The imaging area of the DR3, from which Harris et al. derive their sample, is not the 
same as the spectroscopic area in the DR4. Therefore, for the purposes of this comparison, 
we removed all stars not found in the area of sky common to the two data sets from their 
respective samples. This left 5,340 objects classified as white dwarfs by Harris et al. that 
could potentially have been recovered by Eisenstein et al. Of these, 2,572 were assigned 
spectral fibers in DR4, and 2,346 were ultimately confirmed by Eisenstein et al. to be white 
dwarfs. 

Since we wish to restrict our analysis to single (i.e., non-binary) DA white dwarfs, we 
removed all stars classified as DA+M stars in either catalog. Unfortunately, given that the 
Harris catalog contains no further information as to the type of WD, we were unable to 
remove the non DA stars and simply compare what remains with the Eisenstein sample. 
Instead, we compute the completeness for all of the WDs, under the assumption — explored 
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more fully below — that DAs, as the largest component of the WD population, dominate the 
selection function. 

Figure H] shows a comparison of the two samples. The open symbols are the complete 
Harris et al. sample (excluding those, as mentioned above, with Vtan < 30 km/s, those not 
in the region of sky covered by spectroscopy, and the DA+M stars). The gray squares lie 
outside the cuts in color-color space imposed by Eisenstein et al. They may have spectra in 
SDSS, but they were not fit by Eisenstein et al., and therefore will not have made it into our 
sample. The filled green circles are the stars that are in Eisenstein et al. In other words, if 
the SDSS spectral coverage of WDs were complete, and Eisenstein et al. recovered every WD 
spectra in SDSS, then all of the open circles would be fill ed. The inside of th e blue box is 



the exclusion region for SDSS's QSO targeting algorithm ([Richards et al.ll2002l ). specifically 
implemented to eliminate WDs from their sample. Note that our sample is more complete 
for the stars outside this region. 

Figure [5] shows the discovery probability as a function of g — i color and g magnitude. 
Darker areas mean a higher probability of discovery, with black indicating that all the WDs 
in the Harris et al. sample in that area of color- magnitude space made it into our sample. 
We have performed a box smoothing to eliminate small scale fluctuations. 

There is a drop off in discovery probability for stars bluer than (7 — z ~ —0.2 at all ap- 
parent magnitudes. This corresponds to the red edge of the exclusion region of the QSO 
targeting algorithm, as noted above. The QSO algorithm is also itself a function of apparent 
magnitude, which accounts for the general decrease at fainter magnitudes in the red half 
of the diagram, and the much steeper drop off between 5? ~ 19 and g ^ 19.5. The bluer 
stars {g-i< -0.2), most of which are targeted by the HOT_STANDARD or SERENDIP- 
ITY_BLUE algorithms, show the opposite: a slight increase at fainter magnitudes. 

To give a better sense of the order of magnitude of our completeness. Figure |6] shows a 
histogram of the values in Figure [51 For most of the cells that end up in the bins for 0, 1, 
and 0.5, the Harris et al. sample contains only one or two stars. The mean completeness for 
the whole sample is ~ 51%. 

To derive our final completeness correction, we must further consider the incompleteness 
and contamination in the Harris et al. sample itself. Assuming that the SDSS photometric 
database is essentially complete down to g = 19.5, then the i ncompletenes s in H arris et 



al. comes mainly from two sources: 1) the incompleteness in the [Munn et al.l ([20041 ) proper 
motion catalog, and 2) the tangential velocity limit of 30 km/s imposed, which results in 
some low tangential velocity WDs being dropped from the sample. However, with one 
negligible exception, none of the criteria used to target objects for spectroscopy in SDSS, 
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nor those used by Eisenstein et al. to select white dwarf candidates, depends exphcitly on 
proper motion or tangential velocity. Thus we assume that the low-velocity stars — dropped 
from the Harris et al. sample — will be recovered by Eisenstein with the same probability as 
the high- velocity stars — i.e., the stars in Figure HI 

Contamination poses a bit more challenging problem. At first glance, it would seem 
that the reverse of the above process could be applied, whereby those objects in Harris et 
al. which did get spectral fibers — but were ultimately rejected as WDs by Eisenstein et al. — 
could be removed from the sample, and those that did not get spectra could be assumed to 
follow the same distribution. This latter assumption, however, is unlikely to be true. SDSS 
gives very low priority to targeting white dwarfs specifically, and we would thus expect a 
larger fraction of the objects that get spectral fibers to turn out to be contaminating objects 
(in particular QSOs, of which we found 13 in the Harris et al. sample) than if the fibers 
were assigned purely randomly. Furthermore, many of the 225 objects which have spectra 
in DR4 but are not included in the Eisenstein catalog may actually be white dwarfs which 
Eisenstein's algorithms dropped for some other reason, e.g. they lie outside the color and 
magnitude ranges used for initial candidate selection, or there is a problem (low S/N, bad 
pixels) with the spectrum. Approximately 100 appear to be DC white dwarfs to which the 
SDSS spectroscopic pipeline assigned erroneous redhifts on the basis of weak noise features. 
Ultimately, we have chosen to adopt the contamination fraction of Harris et al. (2%) for the 
whole sample, and have reduced our final completeness correction accordingly. This choice 
has a negligible effect on the small scale structure of the WDLF in which we are interested. 

Finally, we note that the Harris et al. sample has an apparent magnitude limit of 
g = 19.5, whereas the spectroscopic sample contains objects down to g ^ 20.5. Given that 
the SDSS targeting algorithms are themselves functions of apparent magnitude, our com- 
pleteness correction is as well. An extrapolation of our discovery probability is problematic 
in this area, though, because this is just the apparent magnitude where the QSO targeting 
algorithm drops off rapidly. We have decided to impose a magnitude cutoff of = 19.5 in our 
sample. This reduces our sample by nearly a half, with a corresponding increase in counting 
error. However, because SDSS spectra have a small range of exposure times (45-60min), 
fainter apparent magnitude usually translates directly into lower S/N and larger errors in 
derived parameters. 

Figure [7] shows the luminosity functions we derive for different choices of limiting mag- 
nitude. We take the generally good agreement between the curves to indicate that our 
completeness correction is doing its job correctly in the g magnitude direction. 

Figure [8] similarly shows the mass functions we derive for different choices of limiting 
magnitude. In the case of the mass function, the S/N of the spectra becomes a much bigger 
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factor. As a consequence of the essentially constant exposure times of SDSS spectra, the 
parameters [T^s and log g) determined from the spectra of fainter objects have larger errors, 
which causes a larger error i n mass. Thus , the M F is broadened when stars with g > 19.0 
are included. For this reason, iKepler et al.l (120071 ) limited their mass functions to stars with 
g < 19.0, and we follow their lead for the remaining MFs in the current paper. 



5. Luminosity Functions And Discussion 

Figure [9] shows the WD mass function we derive for all stars with T^s > 12, OOOK and 
g < 19.0. The red dashed line is the MF corrected only by 1/Knax — i-e., before we apply 
our completeness c orrection. It generally shows good agreement with the MF derived in 



Kepler et al.l (120071 ) (blue points), not surprising considering we use nearly the same data 
set and very similar WD models. The small differences are due to our use of slightly different 
sets of data and models, as well as differing treatment of duplicate spectra, and can largely 
be considered statistical fluctuations. We refer the interested reader to their paper for a 
more in depth analysis of the WDMF. 

The solid black line in the upper panel shows our MF after correcting for the complete- 
ness of the spectroscopic sample. This curve represents the true local space density of WDs 
per cubic parsec per Mq interval. The bottom panel shows the total weight of each bin 
above the uncorrected MF — essentially the final completeness correction for each bin. There 
is little small scale variation from bin to bin, and our completeness correction mainly has 
the effect of raising the normalization of the whole MF by a factor of ~ 2.2. In other words, 
the shape of the MF is not strongly affected by the completeness correction. 

Figure [10] is th e WDMF for all sta rs down to 8,000K. The dashed red line is for the 



data as reported by iKepler et al.l (120071 ). the dotted blue line is after our correction for the 
upturn in log g. The solid black line is the WDMF for only those stars above 12,000K (i.e., 
the same as Figure [H]) renormalized to the same scale for comparison purposes. There are 
more high mass stars in general, and one spurriously large bin, but on the whole, our log g 
correction recovers a reasonable mass distribution for stars cooler than 12,000K. 

Figure [TT] shows the luminosity function we derive for all of the DA stars in our sample 
down to 7,000K for all stars with g < 19.5. In red is the LF for the data as reported; in 
black is the LF for the data with the increase in log g at low temperature removed. The 
process of removing the excess log g pushes stars to lower masses, making them larger and 
therefore brighter for the same Teg. In the range plotted, the black curve contains a total of 
3,358 WDs, while the red contains 2,940. 
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The lack of agreement between our best LF (black) and the [Harris et al.l (120061 ) lumi- 
nosity function (blue) can be attributed, at least in part, to the differing assumptions used 
in creating the two LFs. Harris et al. derived their temperatures by fitting Bergeron models 
to the photom etry assuming a log q of 8.0 for every s tar, a poor assumption for more than 
30% of WDs (ILiebert et al.l l2005l : iKepler et al.l 120071 ). The temperatures they derive are 
systematically different from the spectroscopic temperatures; Figure [12] shows the fractional 
difference between the spectroscopically and photometrically derived effective temperatures. 
When we use the photometrically derived temperatures and set log g = 8.0, we recover the 
Harris et al. LF fairly well. 

It should also be noted that the Harris et al. luminosity function is for WDs of all types, 
whereas ours is comprised on l y of t he DAs. For each bin in the Harris et al. LF, we have 
used the full lEisenstein et al.l ( l2006l ) catalog to determine a rough DA fraction, and reduced 
the LF reported of Harris et al. accordingly. This DA fraction — s hown in Table 1 — is in 
generally good agreement with previous works (IFleming et al.lll986l ). but we have made no 
attempt to address selection biases in the Eisenstein et al. catalog. 

One other source of the discrepancy between our results and Harris et al. is due to our 
assumption that whatever causes the observed upturn in log g in the cooler stars affects 
only the log g determination and does not alter the spectroscopically derived Tes- As the 
effects of the two parameters on the line profiles are interdependent, this assumption is 
probably not valid. The curves in Figure [TT] suggest that in addition to the excess log 
g, the temperatures determined by line fitting for the cooler stars are probably too high. 
Ultimately, this area of the spectroscopic WDLF will remain uncertain until the problems 
with the model atmospheres have been resolved. 



The LF of iLiebert et al.l (120051 ) shown in green in figure [TT] was compiled from a small 
dataset (348 DA white dwarfs) based on a survey done on photographic plates over 20 years 
ago on a 0.5m telescope. In addition to low number statistics, the dataset suffers from a 
very difiicult-to-quantify incompleteness on the faint end, which is probably responsible for 
the lack of agreement below Mboi ~ 9.5. 



Conclusions 



Our eventual goal is to take advan tage of the tremendous number of WDs spectroscopi- 
cally observed by SDSS and studied by lEisenstein et al.l (120061 ) and others to create separate 
WD luminosity functions for two or more different ranges of mass. This will effectively add 
a third dimension, currently unexplored, to observational WD luminosity functions. 
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In order to carry out this analysis, we must fully understand the manner in which white 
dwarfs w ere selected t o recei ve spectra in SDSS. By comparing the proper -motion selected 



sampl e of iHarris et al.l (I20061) with the spectroscopically determined WDs of iKleinman et al. 



(120041 ) and lEisenstein et al.l (l2006l ). we have derived a WD selection probability over a range 
of parameters that includes nearly the entire useful range of g — i color {—1.0<g — i<0.2) 
and apparent g magnitude (15 < g < 19.5). 

We have also presented additional arguments that the observed upturn in log g is an 
artifact of the model atmosphere line-fitting procedure, or — more likely — a problem with the 
line profiles themselves. Since it may be some time before this problem is fully understood 
and addressed, we have implemented a procedure to remove the excess log g empirically and 
shown that the mass function recovered for the stars cooler than 12,000K reasonably agrees 
with the MF for the hotter stars, which in turn agrees well with previous work. 

Finally, we have presented the first WDLF for spectroscopically determined WDs in the 
Fourth Data Release of the SDSS. In addition to addressing the issues of completeness and 
the observed log g upturn in a more systematic manner than previously attempted, our LF 
contains the largest sample of spectroscopically determined WDs to date (3,358), more than 



than the 298 stars included in the LF of 



six times the 531 presented in lHu et al.l ( 20071). and m o re tha n an order of magnitude more 



Liebert et al.l ((20051). 
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Fig. 1. — log g V. log Teff for the white dwarfs in our sample. At temperatures below 
~12,500K, the log g values begin to rise to an extent unexplained by current theory. The 
solid line is a function empirically fit to the real data. The dashed line is the modest rise 
predicted by theory. The excess at a given Tgfj is subtracted from the measured log g value 
for some of our luminosity functions. 
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Fig. 2. — log g V. log T^s with the upturn in log g removed. 
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Fig. 3. — A comparison of the theoretical colors of the SDSS WDs, derived from the atmo- 
spheric fits (black triangles), with the observed colors, as measured by the SDSS photometry 
(open blue circles). In the upper panel, the colors of the model atmospheres do not agree 
with the observed colors at low temperatures, indicating a problem with the line fitting for 
stars cooler than ~ 12,500K. In the lower panel, where the excess log g has been removed, 
the colors agree much better. 
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Fig. 4. — Color-color plot of the white dwarfs in the two samples used to derive our complete- 
ness correction. Open symbols are WDs from the Harris et al. (2006) sample that a) were 
in the area of sky covered by spectroscopy in DR4, b) had Vtan > 30km/s, and c) were not 
determined by i- and 2;-band excess to be WD -|- n aain-sequence binaries. The filled circles 
are the stars for which SDSS obtained spectra and lEisenstein et al.l (120061 ) confirmed to be 
WDs. The dashed box shows a two-dimensional projection of the QSO targeting algorithm's 
exclusion region. The open gray squares are the WDs from Harris et al. that lie outside 
Eisenstein et al.'s color-color cuts. For clarity, only half of the points have been plotted. 
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Fig. 5. — A map of our completeness correction. Darker areas indicate more complete regions 
of the figure, with black being 100% complete. The overall completeness is of order ~50%. 
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Fig. 6. — A histogram of the completeness values in Figure [51 Most of the 0, 1, and 0.5 
values come from color-magnitude regions in which there are only one or two stars in the 
Harris et al. sample available for comparison. 
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Fig. 7. — Luminosity functions for three different limiting magnitudes. We take the good 
agreement between the curves to indicate that our completeness correction (and the 1/Vmax 
correction) are working properly. 
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Fig. 8. — Mass functions for three different limiting magnitudes. Because of tlie essentially 
fixed integration time for SDSS spectra, objects with fainter apparent magnitudes generally 
have lower signal to noise, which translates directly into larger uncertainties in the derived 
parameters {T^s, log g, and mass). Hence, as we include stars with fainter apparent magni- 
tudes, more stars scatter out of the peak, broadening the mass function. 
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Fig. 9. — The white dwarf mass function for all WDs with Tgfj > 12,000K and g < 19.0. The 
dashed line in the upper panel is the MF corrected only for 1 /Knax, without our completeness 
correction applied. It agrees very well with Kepler et al. (2007 — dots). The solid line is with 
our completeness correction applied, and represents the true local space density of white 
dwarfs. The bottom panel shows the ratio of our two mass functions — i.e., the cumulative 
completeness correction for each bin. The small variation indicates that the completeness 
correction, while changing the overall normalization by roughly a factor of 2.2, has little 
effect on the shape of the MF. 
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Fig. 10. — White dwarf mass functions for WDs with T^g > 8,000K and g < 19.0 both with 
and without the upturn in log g for cooler stars removed. The solid line is the MF from 
Figure [9] renormalized for comparison purposes. 




Fig. 11. — LFs derived in this paper. Removing the log g upturn makes each affected star 
less massi ve, and therefore la rger and brighter, pushin g it to a more leftward Mboi bin. The 
results of Harris et al. J2006h and lLiebert et al.l (120051 ) are shown for comparison. 
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Fig. 12. — A compari son of the sp e ctrall y and photometrically derived temperatures for the 
WDs common to the iHarris et al.l (120061 ) and Eisenstein et al. (2006) samples. 
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Table 1: The fraction of stars in 



Eisenstein et al.l (120061 ) listed as DA or 



DA_auto. Though they generally agree 
with previous results, they should be used 
with much caution, as they were calculated 
crudely and we have taken no care to correct 
for biases in the sample. We have employed 
them here simply to compare our DA-only 
luminosity function to previous work. 



Mm 


DA Fraction 


7.25 


0.9338 


7.75 


0.9243 


8.25 


0.9246 


8.75 


0.8980 


9.25 


0.8433 


9.75 


0.8146 


10.25 


0.7958 


10.75 


0.8158 


11.25 


0.7957 


11.75 


0.7721 


12.25 


0.7985 


12.75 


0.7976 


13.25 


0.8173 


13.75 


0.8009 



